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ABSTRACT 



Models of rotating single stars can successfully account for a wide variety of observed stellar phenomena, such as the surface enhance- 
ments of N and He observed in massive main-sequence stars. However, recent observations have questioned the idea that rotational 
mixing is the main process responsible for the surface enhancements, emphasizing the need for a strong and conclusive test for 
rotational mixing. 

We investigate the consequences of rotational mixing for massive main-sequence stars in short-period binaries. In these systems the 
tides are thought to spin up the stars to rapid rotation, synchronous with their orbital revolution. We use a state-of-the-art stellar 
evolution code including the effect of rotational mixing, tides, and magnetic fields. We adopt a rotational mixing efficiency that has 
been calibrated against observations of rotating stars under the assumption that rotational mixing is the main process responsible for 
the observed surface abundances. 

We find that the primaries of massive close binaries (M\ ~ 20 M Q , P ,b £j 3 days) are expected to show significant enhancements in 
nitrogen (up to 0.6 dex in the Small Magellanic Cloud) for a significant fraction of their core hydrogen-burning lifetime. We propose 
using such systems to test the concept of rotational mixing. As these short-period binaries often show eclipses, their parameters can 
be determined with high accuracy. 

For the primary stars of more massive and very close systems {M\ ~ 50 M Q , P OTb < 2 days) we find that centrally produced helium is 
efficiently mixed throughout the envelope. The star remains blue and compact during the main sequence evolution and stays within 
its Roche lobe. It is the less massive star, in which the effects of rotational mixing are less pronounced, which fills its Roche lobe 
first, contrary to what standard binary evolution theory predicts. The primaries will appear as "Wolf-Rayet stars in disguise": core 
hydrogen-burning stars with strongly enhanced He and N at the surface. We propose that this evolution path provides an alternative 
channel for the formation of tight Wolf-Rayet binaries with a main-sequence companion and might explain massive black hole binaries 
such as the intriguing system M33 X-7. 

Key words, binaries: close - Stars: rotation - Stars: abundances - Magellanic Clouds - Stars: Wolf-Rayet - X-rays: binaries 
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The rotation rate is considered as one of the main initial stel- 
lar parameters, along with mass and metallicity, which deter- 
mine the fate of single stars. Rotation deforms the star to an 
oblate shape (of which Acherna r is an extreme example, see 



(lMaederlll987t lYoon & Langerll2005l) . Some over-luminous O 
and in the Magellanic Clouds in the Magellanic Clouds 

have been put forward as examples pf this type of evolution 



(Bouret et al 



Marti ns et al 



2003 h IWalborn et alJ l2004t iMokiem et afl l2007b 



2008). This alternative evolutionary scenario has 



iDomiciano de Souza et al.l 2003|). it interplays with the mass 
loss from the star (e.g. iFriend & Abbottll 19861 lLangerl Il998t 
iMaeder & Mevnefll2000aF and it induces instabilities in the in- 
terior leading to turbulen t mixing in otherwise stable layers (e.g. 
IMeynet & Maederl Il997l) . As rotationally induced mixing can 
bring processed material from the core to the surface, it has 
been proposed as explanation for observed surface abundance 
anomalies, such as a nitrogen enrichment found in several mas- 
sive m a in-sequence stars (e.g. | Walbornll976tlMaeder & MevneH 
l2000bt iHeger & Langerl l2000h . In models of rapidly rotating 
massive stars, rotational mixing can efficiently mix the centrally 
produced helium throughout the stellar envelope. Instead of ex- 
panding during core H burning as non-rotating models do, they 
stay compact, become more luminous and move blue-wards in 
the Hertzsprung-Russell diagram. This type of evolution is com- 
monly referred to as (quasi-)chemically homogeneous evolution 



been proposed as a way to create rapidly rotating massive helium 
stars as the possible pr ogenitors of long gamma-ray bursts within 
the collapsar scenario (lYoon & Langerll2005UWoosiey & Hegerl 
l2006l:ICantiello et alJl2007l ). 

Multiple attempts hav e been made to constrain the efficiency 
of rotational mixing (e.g. iGies & Lamberllll992t iFliegner et all 



ng 

1996t iDaflon et all l2Q0lt I Venn et al.l l2002h iKorn et al.l I2002: 



Huang & Giesl l2006b iMendel et all l2006l) . but often these at- 



tempts remained inconclusive due to limited sample sizes or a 
strong bias towards stars with low projected velocities. The VLT- 
FLAMES survey provided, for the first time, a large sample of 
massive stars with accurate abundance determinat ions, cover- 
ing a wide range of proj ected rotational velocities dEvans et al.l 
l2005UHunteretalJl2008l) . This sample was used to calibrate the 



1 Wolf-Rayet star s with evidence of e nhanced nitrogen and hydrogen 
in their spectra (e.g. Schnurr et al. 2003) 
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efficiency of rotational mixing, under the assumption that rota- 
tional mixing i s the main proce ss responsible for the observed 
enhancements (IBrott et al.ll2QQ9t) . 

These authors performed a population synthesis, based on 
detailed models of rotating single stars, to reproduce the prop- 
erties of the VLT-FLAMES sample. They find that their models 
cannot account for a large number (20% of the sample) of highly 
nitrogen-enriched, slow rotators and a group of evolved fast ro- 
tators which are relatively non-enriched (a further 20% of the 
sample). This raises the question whether other processes play 
an important role in explaining the nitrogen enhancements of 
massive main-sequ ence stars, such as mass transfer in binaries 
dLanger et al.ll2008l) . 

Regardless of the successes of r otating stellar models to ex - 
plain a variety of stellar phenomena dMaeder & Mevnetll2000bl) . 
rotational mixing is still a matter of debate. Clearly a conclusive 
observational test for the concept of rotational mixing is needed. 
In this paper we propose to use eclipsing binaries for this pur- 
pose. 

Eclipsing binaries have frequently been used to test stellar 
evolution models as they provide accurate stellar masses, radii 
and effective temperatures. Even beyond our own Galaxy, in the 
Magellanic Clouds, stellar parameters of O and early B stars 
have been determined wi th accuracies of 10% (lHarries et alj 
120031: iHilditch et al.ll2005l). which have bee n used to test binary 
evolution models (e.g. lDe Mink et ai1l2007l) . As rot ational mix- 
ing is more important in more massive stars (e.g. IHeger et all 
2000), it is a major advantage to know the stellar masses for 
quantitative testing of the efficiency of rotational mixing. 

In close binaries with orbital periods P rbit less than a few 
days, the tides are so strong that the stars rotate synchronously 
with the orbital motion: P sp j n = P or bit- With stellar radii known 
from eclipse measurements, this enables us to determine the ro- 
tational velocity directly from the orbital period. This is the sec- 
ond important advantage of using binaries for testing rotational 
mixing with respect to single stars. For single stars fitting of 
spectral lines allows only for the determination of v sin i, where 
u is the rotational velocity at the equator and i the inclination of 
the rotation axis, which is generally not known. 

Here, we propose to use eclipsing binaries, consisting of two 
detached main-sequence stars. Detailed calculations of binary 
evolution show that if one of the stars fills its Roche lobe dur- 
ing the main sequence, it does not detach again before hydrogen 
is exhaust ed in the core, except maybe for a very sho rt thermal 
timescale dWellstein et al.l 120011: IDe Mink et al.ll2007l) . Turning 
this around we find that, in a binary with two detached main- 
sequence stars, we can safely exclude the occurrence of previ- 
ous mass transfer. In other words, the fact that a system consists 
of two detached main-sequence stars, constrains the evolution- 
ary history. In contrast, for a fast rotating apparently single star 
we do not know whether the star was born as a rapidly rotating 
single star, or whether the rapid rotation is the result of mass 
transfer in a binary or of a binary merger. The companion, if still 
present, may be very hard to detect, being a faint low-mass star 
in a wide orbiQ. 

If the spectra of a binary are of high quality, one can de- 
termine the surface abundances of the two components (e.g . 
lLeushinlll988tlPavlovski & Hensbergel2005t iRauw et alj|2005l) . 
These surface abundances, together with accurate determina- 
tions of the stellar parameters and the orbital period have the po- 



2 Abundance determinations of boron, if available, may be used 
to distinguish betwe en binary effects and pure rotational mixing, see 
iFliegner et all d 19961) 



tential of strongly constraining the efficiency of rotational mix- 
ing. 



The evolution of close massive binaries has been mod- 



eled by various groups (e.g. iPodsiadlowski et al.l 


1992; 


iPolsl 19941 Wellstein et all 200 it IWellstein & Langeri 


1999; 


Nelson & Eggleton 2001; Belczvnski et al. 20021 Petro vie et al. 



l2005at IVanbeveren et al.H2007l: IDe Mink et al.ll2007l and refer- 



ences therein). In this work we use a detailed binary evolution 
code to predict the surface abundances for massive detached 
close binaries. In addition we discuss models of very massive 
close binaries, in which rotational mixing can be so efficient that 
the change in chemical profile leads to changes in the stellar 
structure. 

2. Stellar evolution code 

We model the evolution of rotating massive stars u sing the ID 
hydrod yna mic stellar evolution co de described by lYoon et all 
(120061) and IPetrovic et ail d2005bl) . which includes the effects 
of rotation on the stellar structure and the transport of angu- 
lar momentum and chemic al species via rota tionally induced 
hydrodynamic instabilities dHeger et al.l l2000b . The rotational 
instabilities considered are: dynamical shear, secular shear, 
Eddington- Sweet circulation a nd the Goldreich-Schubert-Fricke 
instability (IHeger et alj |2000|) . Two processes dominate rota- 
tional mixing in massive stars: Eddington-Sweet circulation, 
large-scale meridional currents resulting from the thermal im- 
balance between pole and equator characteristic of ro tating stars 
(Ivon Zeipellll924t lEddingtonll 19251 [19261: IVogt|[l925l) and shear 
mixing, eddies, that can form between t wo layers of t he star ro- 
tating at different angular velocity (e.g. IZahdfl974l) . We take 
into account an gular momen tum transport by magnetic torques 
as proposed by Spruitl d2002l) as these can successfully provide 
the coupling between the stellar core and envelope necessary to 
explain the observed spins in young compact stellar rem nants 
(IHeger et al.ll2005l: IPetrovic et al1l2005bl: ISuiis et al.|[2Q08h . 

Mixing of chemical species and the transport of angular mo- 
mentum are implemented as diffusion processes. For convection 
we assume a mixing length parameter omlt = 1-5. In semi- 
convective regions we a ssume efficient mixing (asEM = 1 -0, as 
defined in lLangerlll991l) . The turbulent viscosity is determined 
as the sum of the convective and semi-convective diffusion coef- 
ficients and those that arise from the rotationally induced insta- 
bilities. The inhibiting effect of gradients in the mean molecular 
weight on rotatio nal mixing is decreased by a factor = 0.1 
(lYoon et al.ll2006b. 

Brott et all (120091) calibrated the efficiency of rotationally in- 
duced mixing using data from the VLT-FLAMES survey of mas- 
sive stars, assuming that rotational mixing is the main process 
responsible for the observed enhancements. To reproduce the 
extension of the main sequence in the Hertzsprung-Russell di- 
agram, they assumed overshooting of 0.355 times the pressure 
scale height. The main parameter responsible for the efficiency 
of rotational mixing is f c , defined as the contribution of rota- 
tionally induced in st abilities to the diff usion of chemical species 
(IHeger et al.ll2000l) . IBrott et all (120091) found that f c = 0.0228 
is needed to reproduce t he spread in N abu ndances observed in 
the VLT-FLAMES data dHunter et al.ll2008l) . We do not consider 
possible mixing by instabilities related to magnetic buoyancy 
due to windi ng up of the magnetic field lines through differen- 
tial rotation ([ Spruit 2002) as this leads to too efficient mixing 
dHunter et al.ll2008l) . 

Metallicity-depende nt mass loss i n the f orm of stellar wind s 
has been included as in lYoon et alj (120061) : IBrott et al.l d2009l) . 
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For the associated angular momentum loss we assume that mass 
is lost with the specific angular momentum equal to the latitudi- 
nally averaged specific angular momentum of the surface layer. 
The mass loss is thus assumed to be independent of latitude. 

The effect of mass and angula r momentum loss on the bi - 
nary orbit is computed according to Podsiadlows ki et all (1 19921) . 
with the s pecific angular momentum o f the wind calculated ac- 
cording to lBrookshaw & Tavanil (|1993|) . We assume that the or- 
bit is circular and that the spins of the stars are perpendicular 
to the orbital plane. T idal interaction is modeled as described in 
iDetmers et al.l (2008) usi ng the timescale for synchronization by 
turbulent viscosity dZahnl 19771) . see also Section[3] Angular mo- 
mentum is transported between the surface layers and the inner 
regions by magnetic torques and rotational instabilities. At the 
onset of central H burning, we assume that the rotation of the 
stars is rigid and synchronized with the orbital motion. 

Initial Composition 

We discuss single stellar models starting with three different 
initial compositions, representing the composition of the Small 
and Large Magellanic Cloud (abbreviated as SMC and LMC) 
and th e Galaxy (GAL) following the approach by iBrott et all 
(120091) to which we refer for details. For C, O, Mg and Si 
we use the average a bundances determi ned for stars in the 
VLT-FLAMES survey dHunter et alj|2007l) . For the Magellanic 
Clouds we adopt the N abundance measured for HII regions 
which is in agreement with t he lowest observed N abundances 
in the VLT-FLAMES sample dHunter et al.l2007l and references 
therein). F or the remaining heav y elements we adopt solar abun- 
dances by lAsplund et al.l d2005l) scaled down by 0.7 dex for the 
SMC and 0.4 dex for the LMC. For our binary models we have 
assumed the SMC comp osition. The opacity tables of OPAL 
(EH esias & Rogersll 19961) are adopted, where the Fe abundance 
is used to interpolate between the tables of different metallicity. 
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Fig. 1. Timescales for tidal synchronization. The four dashed 
lines indicate where the timescales for tidal synchronization by 
turbulent viscosity r synCjtur b and radiative dissipation T syn c,rad for 
zero-age main-sequence stars are equal to (one percent of) the 
main-sequence lifetime r ms of the most massive star of a binary 
system, see Sect. [3] for details. In the region below the gray- 
shaded bands, tides quickly synchronize the stellar rotation with 
the orbit. The resulting equatorial velocity in kms -1 for the pri- 
mary star is indicated with contour levels, assuming the radius at 
the onset of hydrogen burning. We have assumed a metallicity of 
Z=0.004 and a mass ratio M2/M1 = 0.75, to be consistent with 
the models presented in Sect. [4] The average velocity measured 
for apparently single stars in the SMC, (v eq ) = 175kms _I , is 
plotted with a dotted line. Very short orbital periods are excluded 
(hashed region) as the stars fill their Roche lobe at zero-age. 



3. Stellar rotation under the influence of tides 

In a close binary, angular momentum and kinetic energy can be 
exchanged between the two stars and their orbit through tides. 
The system tends to evolve towards a state of minimum me- 
chanical energy due to dissipative processes. This is an equi- 
librium state, where the orbit is circular, the spins of the stars 
are aligned and perpendicular to the orbital plane and the stars 
are in synchronous rotation with the orbital motion, such that 
P&pm = Aii-bit- How quickly this equilibrium state is approached 
depends on the efficiency of the processes responsible for the 
energy dissipation. 

The most efficient form of energy dissipation takes place in 
turbulent regions of the star, such as convective layers, where the 
kinetic energy of the large-scale flow induced by the tides cas- 
cades do wn to smaller and sm aller scales, until it is dissipated 
into heat. Zahn] (I1977L 1 19891) estimated the timescale for syn- 
chronization due to this "turbulent viscosity" as a function of the 
ratio q of the mass of the companion star to the mass of the star 
under consideration and on the ratio of the stellar radius R over 
the separation a between the two stars: 

Tsync.turb = /turb <7~ 2 i~\ year. (1) 

The proportionality factor / tul b depends on the structure of 
the star, on the location of the turbulent layers and on the 
timescale for dissipation of kinetic energy, which is uncertain. 



Nevertheless, the dependence of the t ime sc ale on R/a is so steep 
that, approximating f tm \, ~ 1 , |Zahnl (1 19771) showed that this ex- 
pression adequately explains the observed orbital period below 
which tides lead to synchronization. 

In the absence of turbulent viscosity, another dissipative 
process is required in order to have efficient tides. If the star 
does not rotate synchronously, it experiences a varying gravita- 
tional potential, which triggers a large range of oscillations in 
the star. These oscillati ons are damp ed near the stellar surface 
by radiative dissipation. |Zahn| (1 19751) derived the corresponding 
timescale for synchronization: 

//?\-17/2 

Tsync.md = /rad <? ^ ( 1 + tf) 5/6 [~J year, (2) 

Here, M denotes the mass of the star under consideration, I its 
moment of inertia and E% the tidal coefficient, which is sensi- 
tive to the structure of the star, i n particular to the size of the 
convective core: £2 <x C^core/^) 8 dZahnl 1977b . 

The early-type massive stars considered in this work have 
radiative envelopes (except for convective zones just below the 
surface which contain almost no mass, e.g. lCantiello et al1l2009l) 
and synchronization by radiative damping has been proposed to 
be the most efficient dissipation mechanism. However, in fast 
rotating stars turbulence is induced by rotational instabilities 
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and o ne may argue that th e first times cale applies also to these 
stars dToledano et alj2007l) . In addition lwitte & Savoniield 19991) 
show that in some cases resonance locking can contribute to even 
more efficient tidal dissipation. 

In Figure [TJ we compare both synchronization timescales 
with the main-sequence lifetime of the more massive star in a 
binary. For this plot we have used the radii of non-rotating zero- 
age main sequence stars assuming a metallicity of Z = 0.004 
and a mass ratio of M2/M1 = 0.75. For E2 we use an analytic 
approximation bv lHurlev et al.l d2002l). for the Roc he-lobe radius 
we use the analytic expression by Eggletonl d 1983b . The diagram 
shows that in binaries with orbital periods shorter than approx- 
imately 3-5 days, the timescale for synchronization is less than 
one percent of the main-sequence lifetime tms, regardless of the 
actual process responsible for the synchronization. The resulting 
equatorial velocity of the primary star, assuming synchronous 
rotation, is given in kms~' along contour lines. 

For comparison: the average equatorial velocity for appar- 
ently single stars in the VLT-FLAMES survey with masses be- 
tween 7 and 2 5 M is 150 kins' 1 for the LMC and 175 kms _I 
for the SMC dHunter et al.1 120081) . In most binaries tides will 
slow down the rotation of the star, but in the tightest binaries 
the rotation rate is higher than that for average single stars. Note 
that the binary systems of interest in this work, for which the 
components have high enough equatorial velocities to show sig- 
nificant surface abundance changes by rotational mixing (>100- 
150 kms -1 , see Sect. @), are in the region where tides are very 
efficient, according to both prescriptions for the synchronization 
timescale. 

4. Rotational mixing in single star models 

Any element produced or destroyed in the hot interior of the star 
can in principle be used as a tracer of rotational mixing. The 
amount by which the surface abundance of a particular element 
changes depends on the location in the star where the element is 
processed. 

Helium is the most important element synthesized during 
core hydrogen burning. By the time a substantial amount of he- 
lium has been formed in the stellar core, a gradient in mean 
molecular weight has been established at the interface between 
the core and the envelope, which has an inhibiting effect on rota- 
tionally induced mixing. Models of moderately fast rotating stars 
(with an equatorial velocity of about 170 kms -1 and a mass of 
about 20 M Q ) show that helium will only appear at the surface 
towards the very end of the main-sequence evolution. In more 
massive stars, the helium surface abundance may be significantly 
enhanced but this is partly due to the strong stellar winds. This 
makes helium not very suitable for testing the effect of rotational 
mixing. 

A better tracer for rotational mixing during the main- 
sequence evolution is nitrogen, which is produced in the hot in- 
terior layers on a very short time scale when carbon is converted 
into nitrogen by CN-cycling. On a longer timescale (about 1-2 
Myr in the center of a 20 M Q star) the CNO cycle comes into 
equilibrium, leading to additional N at the expense of both car- 
bon and oxygen. Rotational mixing can bring nitrogen to the sur- 
face, resulting in a gradual increase of the nitrogen surface abun- 
dance over the main-sequence lifetime (Fig. Nitrogen can 
be measured in rotating B-type stars in the Magel lanic Clouds 
with accuracies of 0.2-0.25 dex dHunter et al.ll2008l) . Carbon de- 
creases accordingly, but cannot be measured as accurately and is 
therefore less suitable as a tracer of rotational mixing. Other ele- 
ments such as boron can be used, but these are less abundant and 



may be hard to detect, especially in metal-poor environments 
such as the Magellanic Clouds. 

Rotational mixing becomes more efficient in higher-mass 
stars since radiation pressure becomes more important, which 
helps to overcome the entropy barrier at the interface between 
the core and the envelope. Also, the ratio of the timescale for 
meridional circulation with respect to the main-sequ ence life- 
time decreases with increasing mass (lYoon et al.ll2006l) . The dis- 
advantage of using very massive stars as test cases for rotational 
mixing is the uncertainty in their mass-loss rates. As mass loss 
exposes deeper layers of the star, enriched in N and He, it has 
qualitatively the same effect on the surface composition as ro- 
tational mixing. At lower metallicity mass loss in the form of 
a radiatively driven stellar wind is reduced. Therefore, also the 
uncertainty in the mass-loss rate is less important. 

Figure|2]shows the nitrogen surface abundance as a function 
of time for rotating single stellar models of 20 M with three 
different initial compositions, representative for the relatively 
metal-poor composition of the Small and Large Magellanic 
Cloud, and for a mixture representing stars in the Galaxy. The 
equilibrium abundances for N in the center for CN- and CNO- 
cycling are plotted as horizontal lines. The Galactic model with 
an equatorial velocity of 230 kms -1 shows nitrogen enhance- 
ments up to about 0.3 dex, while in the SMC model the N 
abundance increases by up to approximately 0.8 dex. The high 
C/N ratio in the Magellanic clouds is partly responsible for 
this effe ct as it leads to a strong increase of N during CN- 
cvcling dBrott et al 112009b . Indeed, in the VLT-FLAMES survey, 
lHunteretal.ld2008l) find a larger spread in N abundances for the 
stars in the Magellanic Clouds than for ones in the Galactic stars. 
We conclude that the Magellanic Clouds are the most promising 
location for testing rotational mixing, as the effect of rotational 
mixing is most pronounced in the N surface abundances. 



5. Binary models 

In this section we present binary evolution models calculated 
with the same set of input physics as the single stellar models 
presented in Section |4] taking into account the effects of mass 
and angular momentum loss on the orbit and spin-orbit coupling 
by tides. We assume a composition representative of the Small 
Magellanic Cloud, which is relatively metal-poor and has a high 
carbon to nitrogen ratio. We discuss two specific model sets: sys- 
tems with a massive primary, Mi = 20 M , in Sect. 15. II and sys- 
tems with a very massive primary, Mi = 50 M Q , in Sect. 



5. 1 . Massive binaries (20 M + 1 5 M G ) 

With the Small Magellanic Cloud sa mpl e of double-lined ecli ps- 
ing binaries bv lHarries et al.l d2003l) and lHilditch et alj d2005l) in 
mind, which contains 21 detached system^ with orbital peri- 
ods ranging from 1-4 days and primary masses ranging from 
7-23 M , we chose to model the following binary systems. For 
the mass of the primary component we adopt 20 M , for the 
secondary component 15 M . The evolution of both stars is fol- 
lowed starting at the onset of central hydrogen burning, at t = 0, 
until the primary star fills its Roche lobe, at t — ?rl. We adopt 



3 Possibly only 20 systems are detached. For two of the systems an 
alternative semi-detached solution exists. For one of these systems a 
comparison to binary evolution models including the effects of mass 
transfer showed that the s emi-detached solutio n was more consistent 
than the detached solution dDe Mi nk et al. 2 0071) . 
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Fig. 2. Surface nitrogen abundance versus time for single stellar models of 20 M with different (approximate) initial equatorial 
rotational velocities starting with the initial composition of the SMC (left panel), LMC (center panel) and GAL (right panel). The 
nitrogen abundance is given in the conventional units: the logarithm of the number fraction of nitrogen «n with respect to hydrogen 
«h, where the hydrogen abundance is fixed at 10 12 , i.e. log e(N) = log 10 [«N/«H] + 12. 

Table 1. Key properties of the massive binary evolution models (20 M +15 M ) as defined and discussed in Sect. 15. II 
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b The surface mass fraction at the onset of Roche-lobe overflow divided by the initial mass fraction of resp. boron ( 10 B + U B), carbon ( 12 C) and 
nitrogen ( 14 N). 

b The equatorial rotational velocity for the primary star, averaged over time from the onset of hydrogen burning until the start of mass transfer 
(one before last column) and at the onset of Roche-lobe overflow (last column). 



orbital periods up to 3 days. In these systems the tides are effi- 
cient enough to keep both stars in synchronous rotation with the 
orbit (Sect.©. 



Surface abundances 

The abundances of nitrogen, carbon, helium and boron at the sur- 
face of the primary star are shown as a function of time in Fig. [3] 
see also Tab. Q] The nitrogen abundances at the surface starts to 
increase after about 1 Myr for the 1 . 1 day binary, and after about 
5 Myr for the 3.0 day binary. This time delay fdeky is the time it 
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Fig.3. Surface abundances of nitrogen ( 14 N), carbon ( 12 C), boron ( 10 B+ n B) and the mass fraction of helium at the surface versus 
time for a 20 M star with a 15 M Q close companion. Note the different vertical scales. The abundance of an element X is given 
relative to hydrogen in the conventional units: log e(X) = log 10 (nx/«H) + 12, where rix and n# refer to the number fractions. The 
different lines show the evolution assuming initial orbital periods of 1.1 days and between 1.2 and 3 days, with a spacing of 0.2 
days. The tracks are plotted from the onset of central H burning until the onset of Roche-lobe overflow. See also TableQ] 



takes to transport the nitrogen from the deeper layers, where it 
is produced, to the surface^]. The largest enhancement, 0.6 dex, 
is achieved in the 1 . 1 day binary. The shorter the orbital period, 
the faster the rotation of the stars, the more efficient rotational 
mixing and the faster the surface abundances change with time. 
On the other hand in the systems with short orbital periods the 
stars fill their Roche lobe at an earlier stage, leaving less time to 
modify their surface abundances. 

The typical uncertainty in the nitrogen abundance determi- 
nations for stars in the VLT- FLAMES survey is about 0.2 dex. 
Therefore, to be able to detect whether nitrogen is enhanced, it 
should exceed the initial nitrogen abundance by 0.2 dex. The 
age at which this lower limit for detecting a surface nitrogen en- 
hancement is exceeded is denoted as t el ± in Table Q] Systems 
with orbital periods shorter than 3 days all reach surface en- 
hancements of 0.2 dex, at the latest just before they fill their 
Roche lobe. The time span during which the primary can be 
observed with a surface nitrogen abundance exceeding 6.7 dex, 
listed as fRL-fenh in TableQ] is highest for the 1.4 day system, 1.5 
Myr. For detached tidally locked main-sequence binaries with 
initial orbital periods greater than about 3 days we do not expect 
detectable N enhancements on the basis of these models: their 



4 We define r^iay as the age at which the surface nitrogen abundance 
is enhanced by 0.01 dex. 



rotation rates and therefore the efficiency of rotational mixing is 
too low. 

While the nitrogen abundance increases, the carbon abun- 
dance decreases accordingly, see Fig. [3] It decreases by less than 
0.1 dex in our models. The changes in the mass fraction of he- 
lium at the surface are very small, just over 1% at maximum 
for the tightest system. An element that is a very sensitive tracer 
of rotational mixing is boron. It is easily destroyed in the hotter 
layers just below the surface. Therefore the time delay for boron, 
after which the boron surface abundance starts to change, is very 
short. The change in the boron surface abundance is consider- 
able, up to 1.2 dex in the 1.4 day model. In practice it may be 
hard to measure boron due to its low overall abundance. 

The closer the stars are to filling their Roche lobe, the higher 
is their surface N abundance. This can be turned around and may 
be used to predict which binary systems are likely to show nitro- 
gen surface enhancements. In Tab. Q] we indicate the Roche-lobe 
filling factor R/Rrl, defined as the radius of the primary star 
over the Roche-lobe radius, at the moment when the surface N 
enhancement exceeds the initial abundance by 0.2 dex. Based on 
these models, we predict that the observed systems with masses 
close to 20 and 15 M have to fill their Roche lobes by about 
86% or more before they show detectable surface N enhance- 
ments (see TableQ]col. 5). We do note that the radii predicted by 
our models are sensitive to the assumed amount of overshooting. 
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In the last two columns of TableQ]we indicate the rotational 
velocity at the equator, averaged over time and at the onset of 
Roche-lobe overflow. Whereas a single star at this metallicity 
slows down its rotation rate due to its evolutionary expansion, 
the spin periods of the stars in tidally locked binaries are nearly 
constant because the orbit acts as an angular momentum reser- 
voir. Therefore rotationally induced mixing is more important in 
a star in a tidally locked binary than in a single star that started 
with the same initial rotation rate. The effect is small, however, 
as the expansion mainly takes place after about 5-6 Myr for a 
20 M star. Binaries with orbital periods shorter than 2 days 
fill their Roche lobe before the expansion sets in. The effect is 
stronger at higher metallicity, where angular momentum loss by 
the stellar wind becomes important, see Sect. [6] 

Expected trends with system mass and mass ratio 

For an ideal comparison between models and observations one 
would prefer a model in which the binary parameters match the 
observed parameters. In the above, we discussed models with a 
specific choice for the primary mass and mass ratio and we only 
varied the initial orbital period. Here we discuss how our find- 
ings are expected to change for systems with slightly different 
system masses and mass ratios. 

In more massive binaries rotational mixing is more efficient, 
such that they will show a shorter time delay, even when mea- 
sured in units of the core hydrogen-burning lifetime. Detectable 
surface enhancements are expected to occur at an earlier stage, 
while the stars are further away from filling their Roche lobe. 
At very high masses additional effects can play a role, which 
are discussed in Sect. 15.21 The secondary star in the models 
discussed above does not show significant surface abundance 
changes for C, N and He. The surface boron abundance does 
decrease by up to 0.6 dex for the 1 .4 day model. In systems with 
mass ratios closer to 1, the time scale for rotational mixing be- 
comes similar in both stars and nitrogen enhancements are ex- 
pected for both stars. 

Another effect of changing the mass ratio is that the size of 
the Roche lobe changes. A more extreme mass ratio results in 
a wider Roche lobe for the primary, if the primary mass and or- 
bital period are kept constant. The efficiency of rotational mixing 
does not change, as it depends on the rotational period which is 
fixed by the orbital period. However, the star has more space 
to expand before it fills its Roche lobe. We noted above that 
the primary stars fill their Roche lobes by at least 86% before 
they show detectable surface N enhancements (see also Sect[6j). 
This changes to approximately 83% and 76% for mass ratios of 
q = 0.5 and q = 0.25 respectively. A bigger Roche lobe also 
leaves more time to enhance the surface abundance. This im- 
plies a greater chance of catching the system in a stage where 
the surface abundance is significantly enhanced. 

To summarize: the biggest surface N enhancements are ex- 
pected for systems with a high-mass primary, which is close to 
filling its Roche lobe, with a companion which has a signifi- 
cantly lower mass in an orbit of no more than a few days. 

5.2. Very massive binaries (50 M +25 M Q J 

In more massive binaries, rotational mixing can be so efficient 
that the change in chemical profile leads to significant structural 
changes. In this model set we adopt a primary mass of 50 M , a 
secondary mass of 25 M and orbital periods varying from 1.5 
to 4 days, assuming an SMC composition to be consistent with 




Fig. 4. The evolution from the onset of central H burning until 
the moment of Roche-lobe overflow for a 50 M star in a binary 
with a 25 M companion (not plotted) with initial orbital periods 
between 1.5 and 4 days. 
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Fig. 5. Helium abundance at the surface as a function of the 
helium abundance in the center for the same systems as plotted 
in Figure |U 



the models presented in Sect. I5.ll Although such massive close 
systems are rare, observational counterparts do ex ist, for exam- 
ple tw o of the four massive binaries presented by iMassev et alj 
(I2002h which are located in the R136 cluster at the center of 
the 30 Doradus nebula in the Large Magellanic Cloud: R136- 
380 and R 136-4^. Another example with an even closer orbit is 
[L72] LH 54-4253 locate d in the LH 54 OB association in the 
Large Magellanic Cloud dWilliams et al.ll2008l) . All three binary 
systems have O-type main-sequence components, which reside 
well within their Roche lobes. 

Figure [4] shows the evolution of the primary stars in the 
Hertzsprung-Russell diagram, until one of the stars in the binary 
fills its Roche lobe (not necessarily the primary, see below). At 
the onset of hydrogen burning their location in the diagram is 
very similar, although the stars in tighter binaries, which rotate 
faster, are slightly cooler and bigger: a direct consequence of 



23.4: 
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Mi = 56.9 ± 0.6 Mq, M 2 
M, = 40.3 ± 0.1 M G , M 2 = 32.6 ± 0.1 M and P mb = 2.89d 
Mi = 47 ± 2 Ms, M 2 = 28 ± 1 M Q and P orb = 2.25d 
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Fig. 6. Nitrogen abundance as a function of age at the surface for the primary and secondary star of the same systems as plotted in 
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Fig. 7. Radius as fraction of the Roche-lobe radius for the primary (left panel) and secondary star (right panel) of the same systems 
as plotted in Figure [4] 



Table 2. Key properties of very massive binaries (50 M +25 M ) as described in Sect. I5.2l (see also Table[T]and Sect. I5.lt . In the 
last column we indicate which component fills its Roche lobe first. 
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the centrifugal force. As they evolve their tracks start to deviate. 
The wider systems (f or b > 2.0d) evolve similarly to non-rotating 
stars: they expand during core hydrogen burning, evolving to- 
wards cooler temperatures until they fill their Roche lobe. Their 
evolutionary tracks overlap in the Hertzsprung-Russell diagram. 

The primaries in the tighter systems (P rb < 2.0d) be- 
have very differently: they evolve left- and up-ward in the HR- 
diagram, becoming hotter and more luminous while they stay 
relatively compact. The transition in the morphology of the 
tracks around P or b * 2 d ays is similar to the bifurcation found f or 
fast rotating sing le stars dMaederHl987llYoon & Langerj|2005l) . 



Surface abundances 

Rotational mixing is so efficient in these systems that even 
a large amount of helium can be transported to the surface. 
Figure [5] depicts the helium mass fraction at the surface as a 
function of the helium mass fraction in the center. In the hypo- 
thetical case that mixing would be extremely efficient throughout 
the whole star, the surface helium abundance would be equal to 
the central helium abundance at all time. This is indicated by 
the dotted line. For the widest systems, the surface helium mass 
fraction is not affected by rotational mixing at all, while for the 
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tighter systems Xn e reaches up to 65%. They follow the evolu- 
tion of chemically homogeneous stars. Figure [5] also shows that 
for each system, mass transfer starts before all hydrogen is con- 
verted into helium in the center. Note that the highest central He 
mass fraction when mass transfer starts is reached in the 1.7 day 
system (more than 80 %), whereas on the basis of stan dard binary 
evolution theory (e.g. lKippenhahn & Weigerj|1967h one would 
expect this to occur in the widest system. This anomalous be- 
havior is connected to the evolution of the radius, as discussed 
below. 

All systems show big enhancements of nitrogen at the sur- 
face of the primary, see Fig. [6] The wider systems are enhanced 
by up to 0.9 dex. In the tight systems the enhancement reaches 
almost 2 dex. This extreme increase is partly due to the fact that 
abundance is measured relative to hydrogen, which is signifi- 
cantly depleted at the surface of the primaries in the tightest 
systems (cf. Table. Also the secondary stars show nitrogen 
surface enhancements, of up to 0.5 dex. 

Evolution of the radius 

The increase of helium in the envelopes of the primary stars in 
the tightest binaries leads to a decrease of the opacity and an in- 
crease in mean molecular weight in the outer layers, resulting in 
more luminous and more compact stars. In Figure|7]we plot the 
stellar radii as a fraction of their Roche-lobe radii. The primary 
stars in the wider systems expand and fill their Roche lobe af- 
ter about 2.5-3 Myr. Contrary to what one might expect, we find 
that Roche-lobe overflow is delayed in tighter binaries. Whereas 
classical binary evolution theory predicts that the primary star is 
the first to fill its Roche lobe, we find instead that for systems 
with P or b < 1.7 days it is the less massive secondary star that 
starts to transfer mass towards the primary. During this phase of 
reverse mass transfer from the less massive to the more massive 
star, the orbit widens. Nevertheless we find that, if we continue 
our calculations, the two stars come into contact shortly after the 
onset of mass transfer and the stars are likely to merge. 

Expected trends with system mass and mass ratio 

In more massive systems the effect of rotational mixing be- 
comes stronger in both components. This may allow for chem- 
ically homogeneous evolution to occur in binary systems with 
wider orbits. If the mass ratio is closer to one, M2 ~ M\, the 
effect of rotational mixing becomes comparable in both stars, 
and they may both evolve along an almost chemically homo- 
geneous evolution track. In that case both stars can stay within 
their Roche lobe and gradually become two compact WR stars 
in a tight orbit. In a system with a more extreme mass ratio, 
M2 <K Mi , the secondary star will hardly evolve or expand dur- 
ing the core hydrogen-burning lifetime of the primary. Also in 
this case Roche-lobe overflow may be avoided during the core 
hydrogen-burning lifetime of the primary component, leading to 
the formation of a Wolf-Rayet star with a main-sequence com- 
panion in a tight orbit. The parameter space in which this type 
of evolution can occur will be examined in a subsequent paper. 

6. Discussion 

We have shown that rotational mixing can have important ef- 
fects in close massive binaries. The mixing parameters in our 
code have been calibrated against observations of rotating mas- 
sive stars in the VLT-FLAMES survey under the assumption that 



rotational mixing is the main process responsible for the ob- 
served surface N enhancements (see Sect. |2|. The predictions 
we present for close binaries can be used to test the validity of 
this assumption. 

In Sect. |4] we conclude that eclipsing binaries in the 
Magellanic Clouds are the most promising test cases for rota- 
tional mixing. In these systems we expect the biggest enhance- 
ments of N at the surface. Also, the uncertainty in the mass- 
loss rates plays a less important role, as radiatively driven winds 
are reduced at low metallicity. A disadvantage of Magellanic 
Cloud binaries with respect to Galactic systems is the greater 
distance. Long integration times may be needed to obtain spec- 
tra with high enough quality to determine the stellar parame- 
ters and abundances. The advantage of using eclipsing binaries 
is that we can compare directly to evolution models with corre- 
sponding masses and orbital periods. Therefore, just a few well- 
studied systems may be enough to put constraints on the effi- 
ciency of rotational mixing. However, ideally one would prefer 
a large sample to enable a statistical comparison. 

The main parameter in our code that affects the efficiency 
of mixing due to rotational instabilities is f c (see SectO which 
has been calibrated directly against the surface N abundances 
in the VLT-FLAMES survey. However, the calibration involves 
multiple parameters, such as which is a measure of how ef- 
fectively a gradient in mean molecular weight can inhibit ro- 
tational mixing. We expect that uncertainties in this parameter 
are not important for our predictions for nitrogen, as this ele- 
ment is mainly produced and transported to the envelope early 
in the evolution, before a strong mean molecular-weight gradi- 
ent has been established at the interface between the core and 
the envelope. However, a lower value of may lead to higher 
helium surface abundances, which may facilitate the possibility 
of chemically homogeneous evolution in close binaries, but this 
remains to be investigated. 

We use a large amount of overshooting in our models (0.355 
times the pressure scale height H p ), to reproduce the exten- 
sion of the main sequence observed in the VLT-FLAMES data. 
However, the amount of overshooting is an uncertain param- 
eter and some authors q uote lower values for the amount of 
overshooting, for example fSchroder et al.l d 1997b who find 0.24- 
0.32H p based on eclipsing binaries with stellar masses between 
2.5 a nd 6.5 M G , (see also IStothers & Chid [19921: lAlongi et al.l 
119931) . We recomputed one of our models (20+15 M , 3 days) 
without overshooting we find that, although the stars are less 
luminous as expected, the surface nitrogen abundance and the 
radius at a given time are very similar (deviations of less 0.01 
dex in the N abundance at a given age, and less than 4% in the 
radius for a given nitrogen surface abundance). For tighter bina- 
ries we expect an even smaller effect on the radius. We conclude 
that our predictions are not very sensitive to the uncertainties in 
the overshooting parameter. 

In our models we assume that mixing processes in binaries 
operate in the same way as in single stars. For both single stars 
and binary members we find that the stellar interior rotates nearly 
rigidly (at least during the early phase of evolution of interest 
here) as a result of efficient internal angular momentum transport 
by magnetic torques. Having very similar internal rotational pro- 
files, the only difference in our models arises from the evolution 
of the rotation rate, which in single stars is governed by angular 
momentum loss and evolutionary expansion, while in binaries 
the tides play a major role. However, in the binary models we 
consider the stars are close to filling their Roche lobe. They are 
slightly deformed in a lob-sided way, no longer being symmet- 
ric around the rotation axis. In addition, one side of the star is 



10 



S. E. de Mink, M. Cantiello, N. Langer, et al.: Rotational mixing in massive binaries: 



irradiated by the companion and may be heated. Although the 
system is synchronized, the tides continue to extract or deposit 
angular momentum from or onto the stars. How such effects, in- 
duced by the presence of the companion star, interact with the 
different rotational instabilities is not well understood and poses 
an additional uncertainty on our predictions. If such effects are 
important and lead to additional mixing, our predictions for the 
N surface abundance can still be used as a test for rotational 
mixing if they are considered as lower limits to the expected N 
abundance. 



Avoiding mass transfer in short-period binaries 

We have shown that rotational mixing, if it is as efficient as as- 
sumed in our models, can lead to chemically homogeneous evo- 
lution for tight binaries with a 50 M primary. In these models 
the primary star stays so compact that the secondary star is the 
first to fill its Roche lobe. 

This peculiar behavior of the radius of stars, which are ef- 
ficiently mixed, has been noted in models of rapidly rotating 
massive single stars (lMaededl!987l) and has been suggested as 
an evo lutionary channel for the progenitors of lo ng gamma-ray 
bursts (lYoon et al.ll2006t IWooslev & Hegerll2006l) in the collap- 



sar scenario 



Wooslevl (I1993F 7 In single stars this type of evo- 



lution only occurs at low metallicity, because at solar metallic - 
ity mass and angular momentum loss in the form of a stellar 
wind spins down the stars and prevents initially ra pidly rotat- 
ing stars from evolving chemically homogeneously dYoon et all 



l2006t iBrott et al.|[2009t) . In a close binary tides can replenish the 
angular momentum, opening the possibility for chemically ho- 
mogeneous evolution in the solar neighborhood. 

The binary models presented here all evolve into contact, but 
(as we discussed briefly in Sect l5.2b Roche-lobe overflow may be 
avoided altogether in systems in which the secondary stays com- 
pact, either because it also evolves chemically homogeneously, 
which may occur if M\ « Mi, or because it evolves on a much 
longer timescale than the primary, when Mi «: Mi. Whereas 
standard binary evolution theory predicts that the shorter the or- 
bital period, the earlier mass transfer sets in, we find that bina- 
ries with the lowest orbital periods may avoid the onset of mass 
transfer altogether. This evolution scenario does not fit in the 
traditional classification of interacting binaries into Case A, B 
and C, based on the evolutionary stage of the primary compo- 
nent at the onset of mass transfer dKippenhahn & Weigert|[l967t 
lLauterbornli970T) . In the remainder of this paper we will refer to 
this new case of binary evolution, in which mass transfer is de- 
layed or avoided altogether as a result of very efficient internal 
mixing, as Case M. 

The massive and tight systems in which Case M can occur 
are rare. Additional mixing processes induced by the presence 
of the companion star, which may be important in such systems, 
will widen the parameter space in which Case M can occur: it 
would lower the minimum mass for the primary star and in- 
crease the orbital period below which this type of evolution oc- 
curs. The massiv e LMC binary IL72] L H 54-425, with an orbital 
period of 2.25 d dWilliams et al.ll2008l see also Sec. 15. 2\ may be 
a candidate for this type of evolution. Another interesting case 
is the galactic binary WR20a, which consists of two core hydro- 
gen burning stars of 82.7 + 5.5 and 81.9 ± 5.5 M Q in an orbit of 
3.69 d. Both stars are so compact that they are detached. The sur- 
face abundance show evidence of rotational mixing: a nitrogen 
abundance of six times solar is observe d and carbon is depleted 
dBonanos et al.ll2004t iRauw et al.ll2005l) . 



Short-period Wolf-Rayet and black-hole binaries 

If Roche-lobe overflow is avoided throughout the core hydrogen- 
burning phase of the primary star, both stars will stay com- 
pact while the primary gradually becomes a helium star and 
can be observed as a Wolf-Rayet star. Initially the Wolf-Rayet 
star will be more massive than its main sequence companion, 
but mass loss due to the strong stellar wind may reverse the 
mass ratio, especially in systems which started with nearly equal 
masses. Examples of observed short-period Wolf-Rayet bina- 
ries with a main-sequence companion are CQ CepQ, CX Ce 
HP 193576E3 and the very massive system HD 31188 _ 
(Ivan der Huchlll200l]) . Such systems are thought to be the result 
of very no n-conservative mass transfer or a common envelope 
phase (e.g. lPetrovic et al.ll2005al) . Case M constitutes an alterna- 
tive formation scenario which does not involve mass transfer. 

Case M is also interestin g in the light of massive black-hole 
binaries. lOrosz et ail d2007l) recently published the stellar pa- 
rameters of M33 X-7, located in the nearby galaxy Messier 33 
which harbors one of the most massive stellar black holes known 
to date, Mbh = 15.7 ± 1.5 M , orbiting a massive O star, Mo = 
70 ± 7 M , which resides inside its Roche lobe in spite of the 
fact that the orbit is very tight, P 01 -b = 3.45 d. The explanation 
for the formation of this system with standard binary evolution- 
ary models involves a common-envelope phase that sets in after 
the end of core helium burning (Case C), as the progenitor of 
the black hole must have had a radius much greater than the 
current orbital separation. This scenario is problematic as it re- 
quires that the black-hole progenitor lost roughly ten times less 
mass before the onset of Roche-lobe overflo w than what is cur - 
rently predicted by stellar evolution models dOrosz et al. 1 120071) . 
An additional problem is that the most likely outcome of the 
common envelope phase would be a merger, as the envelope s 
of massive stars are tightly bound dPodsiadlowski et ai1l2003l) . 
In the Case M scenario the black-hole progenitor can stay com- 
pact and avoid Roche-lobe overflow at least until the end of core 
helium burning, such that it retains its envelope. 

There are examples of stellar mass black-hole binaries with 
short periods and a massive main-sequence companion, in which 
nearly chemically homogeneous evolution may be important. 
IC 10 X-l is a system harboring the most massive stellar mass 
black hole known to date with a mass of at least 21 M , orbiting 
a Wol f-Rayet star of approximately 25 M Q in an orbit of 1 .45 
days dSilverman & Filippenkoll2008l) . Homogeneous evolution 
helps to explain the high mass of the black hole, but the short 
orbital period poses a difficulty, also for Case M: strong mass 
loss during the Wolf-Rayet life time will widen the orbit. Other 
examples of high mass black hole binaries with short orbital pe - 
riods are th e famous systems Cyg X- 10 dHerrero et al.lll995l). 
LMC X- H dOrosz et al.l2008l) and LMC X-30 ( Yao et al.l2005L 
and references therein). 

The subsequent evolution of tight rapidly rotating Wolf- 
Rayet binaries remains to be investigated. If one or both mem- 
bers of the system can reta in enough angu lar momentum to ful- 
fill the collapsar scenario dWooslevl ll993l). which may be hard 
as the tides can slow down the stars (e.g. 



<$), which may be h ard 
Detmers et a"i1l2008l) . it 



CQ Cep: M WR = 24 M Q , M = 30 M , P mb = 1.6 d 
CX Cep: M WR = 20 M Q , M = 28 M G , P mb = 2.1 d 
HD 193576: M WR = 9 M , M = 29 M Q , P orb = 4.2 d 
HD 311884: M WR = 51 M G , M = 60 M , P mb = 6.2 d 

5.6 d 



Cyg X-l: M bh * 10 M G , M * 18 M G , P mb 
LMC X-l: M bh « 10 M G , M Q ~ 30 M G , P B1 



3.9 d 



14 LMC X-3: Mbh * 4-10 M G , M ~ 40 M G , P mb = 4.2 d 



S. E. de Mink, M. Cantiello, N. Langer, et al.: Rotational mixing in massive binaries: 



11 



may lead to the production of one or even two long gamma-ray 
bursts. 

Although promising, the importance of this channel may be 
small due to the limited binary parameter space in which Case M 
can occur. To make any strong statements about this new evo- 
lutionary scenario, further modeling is needed, which we will 
undertake in the near future. 

7. Conclusion 

We investigated the effect of rotational mixing on the evolution 
of detached short-period massive binaries using a state of the 
art stellar evolution code. The efficiency of rotational mixing 
was calibrated under the assumption that rotational mixing is 
the main process responsible for the observed N enhancements 
in rotating stars. 

We find nitrogen surface enhancements of up to 0.6 dex for 
massive binaries in the SMC. The largest enhancements can be 
reached in systems with orbital periods less than about 2 days, in 
which the primary is massive (about 20 M Q or more) and evolved 
(filling its Roche lobe by about 80% or more) and the secondary 
is significantly less massive, which leads to a more spacious 
Roche lobe for the primary (preferably M%jM\ < 0.75). 

We propose to use such systems as test cases for rotational 
mixing. These systems often show eclipses and big radial veloc- 
ity variations, such that their stellar parameters, the rotation rate 
and possibly their surface abundances can be determined with 
high accuracy. This enables a direct comparison between an ob- 
served system and models computed with the appropriate stellar 
and binary parameters. An additional major advantage of using 
detached main-sequence binaries is the constraint on the evolu- 
tionary history. For a fast spinning apparently single star we do 
not know whether it was born as a fast rotator or whether its ro- 
tation rate is the result of mass transfer or merger event. In a de- 
tached main-sequence binary we can exclude the occurrence of 
any mass transfer phase since the onset of core-hydrogen burn- 
ing (see Sect[T|i. 

In the most massive binaries we find that rotational instabil- 
ities can efficiently mix centrally produced helium throughout 
the stellar envelope of the primary. They follow the evolution 
for chemically homogeneous stars: they stay within their Roche 
lobe, being over-luminous and blue compared to normal stars. 
Due to large amount of nitrogen and helium at the surface these 
stars can be observed as Wolf-Rayet stars with hydrogen in their 
spectra. In contrast to standard binary evolution, we find that it 
is the less massive star in these systems that fills its Roche lobe 
first. 

There may be regions in the binary parameter space in which 
Roche-lobe overflow can be avoided completely during the core 
hydrogen-burning phase of the primary. The parameter space for 
this new evolutionary scheme, which we denote Case M to em- 
phasize the important role of mixing, increases if additional mix- 
ing processes play a role in such massive systems. It may provide 
an alternative channel for the formation or tight Wolf-Rayet bi- 
naries with a main-sequence companion, without the need for 
a mass transfer and common envelope phase to bring the stars 
close together. This scenario is also potentially int eresting for 
tight massive black hole binaries, such as M33 X-7 (lOrosz et al.l 
|2007|) . for which no satisfactory evolutionary scenario exists to 
date. 
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